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Abstract. We perform an analysis of preliminary data on strange particles 
yields and fluctuations within the Statistical hadronization model. We begin 
by describing the theoretical disagreements between different statistical models 
currently on the market. We then show how the simultaneous analysis of yields 
and fluctuations can be used to differentiate between the different models, and 
determine if one of them can be connected to underlying physics. We perform a 
study on a RHIC 200 GeV data sample that includes stable particles, resonances, 
and the event- by-event fluctuation of the /C/tt ratio. We show that the equilibrium 
statistical model can not describe the fluctuation, unless an unrealistically small 
volume is assumed. Such small volume then makes it impossible to describe 
the total particle multiplicity. The non-equilibrium model, on the other hand, 
describes both the Kj-n fluctuation and yields acceptably due to the extra boost 
to the TT fluctuation provided by the high pion chemical potential. A(f520) and 
K* abundance is described within error bars, but the S* is under-predicted to ~ 
f .5 standard deviations. We suggest further measurements that have the potential 
to test the non-equilibrium model, as well as gauge the effect of re-interactions 
between hadronization and freeze-out 
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1. Statistical model descriptions of particle yields 

One of the main objectives of heavy ion physics is to study the collective properties 
of strongly interacting matter. It's equation of state, transport coefficients, degree of 
equilibration and phase structure, and the dependence of these on energy and system 
size. 

Thus, the natural approach to study soft particle production in heavy ion 
collisions is through statistical mechanics techniques. Such an approach has a long and 
illustrious history Pl^lElIll- However, the systematic quantitative comparison of data 
to the statistical model is a comparatively recent field; A consensus has developed that 
the statistical hadronization model can indeed fit most or all particles for AGS,SPS 
and RHIC energies □ El COl E] . 

This consensus, is, however, somewhat superficial. The models cited in the 
references above are considerably different w.r.t. each other, beyond their common 
origin in statistical mechanics. As a result, linking fit results to physics in a model- 
independent way has proven to be somewhat problematic. 

For instance, the observation of "the horn", a sharp dependence of certain 
observables, notably the /n'^ ratio, with energy ^21 has generated a lot of 
theoretical interest, and is generally seen as remarkably suggestive of a phase 
transition. However, many completely different explanations were given for this 
phenomenon, all to some extent based on statistical mechanics: 

The original A'+/7r+ "horn" dependence was predicted J3| as a signature of 
deconfinement; The motivating argument was that the onset of the phase transition 
should be marked by a sharp change in the strangeness over entropy ratio, tracked 
by K/tt. This idea has recently been realized [21 in terms of non-equilibrium at 
freeze-out: The extra strangeness and entropy in the high-temperature phase translate 
into lack of chemical equilibrium afler the phase transition, resulting in super-cooling 
(freeze-out temperature ~ 140 MeV) and phase-space over-saturation (over-abundance 
of both light and strange quarks w.r.t. equilibrium expectations). Other authors, 
however, have interpreted the horn without advocating deconfinement; Models have 
been developed where the sharp structure arises due to the transition to the Grand 
Canonical limit for strangeness or from freeze-out conditions where mesons rather 
than baryons carry the bulk of the entropy and conserved quantum numbers |14j . 
Another proposed explanation relies on the absence of equilibration between Kaons 
and pions JHl ■ While this variety of different models is a sign of both the theoretical 
and experimental vitality of our field, it would be very useful to rule some of these 
models out. 

A similar ambiguity is present in the interpretation of strangeness enhancement, 
now firmly observed both at SPS il6^ and RHIC JJTj energies. It was originally 
proposed on the grounds that a perturbative gluon-rich QGP is a lot more efficient as 
a strangeness production mechanism than a hadron gas ^ . However, a considerable 
amount of recent attention has been given to models where the enhancement is 
explained by a suppression of strange particles in low-multiplicity systems due to 
deviation from the Grand-Canonical limit '181. The continued ambiguity over these 
models is evidence that further observables need to be accessed before either of the 
strangeness production mechanisms can be considered "falsified" . 

Finally, a topic of considerable interest is the relationship between the statistical 
model and freeze-out. The RHIC HBT puzzle ^1 1^ makes it clear that freeze- 
out conditions are, at the moment, not well understood. In particular, it is unclear 
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what, if any, is the effect of the evolution between hadronization (the formation of 
hadrons as effective degrees of freedom) and freeze-out (the moment when ah hadrons 
decouple). As shown in |2U| . the discrepancy with experiment would lessen if chemical 
and kinetic freeze-out were to coincide. While, as claimed in [201, such a high freeze-out 
temperature would spoil agreement with particle spectra, fits based on a single freeze- 
out model prescription have shown that, provided a correct treatment of resonances is 
maintained, particle spectra are compatible with simultaneous freeze-out pi, 22, 2^. 
This somewhat invalidates the claims 23| that "blast-wave" fits to spectra (where 
resonance feed-down is not included) are evidence that freeze-out occurs after a long 
period of re-interaction, with S s and f2 s freezing out before tt, K and protons. 

A logical step to clarify these scenarios is to directly measure the abundance of 
resonances. Here, the situation becomes even more ambiguous: As pointed out in 
|25| . resonance abundance generally depends on two quantities: m/T, where m is 
resonance mass and the chemical freeze-out temperature, as well as tF, where V is 
the resonance width and r is the reinteraction time. Observing two ratios where the 
two particles have the same chemical composition, but different m and F, such as 
A(1520)/A vs K*/K, or E(1385)/A vs K*/K, could therefore be used to extract the 
magnitude of the freeze-out temperature and the re-interaction time. Studies of this 
type 123 are still in progress; As we will show, A(1520)/A and K* / K seem to 
be compatible with sudden freeze-out, provided freeze-out happens in a super-cooled 
over-saturated state, as suggested in (SJIOIE]. Other prehminary results, such as /o/tt, 
A/p |H], and now 5]*/A j2Z] seem to be produced in excess of the statistical model 
|2niElj both equilibrium and not. It is difficulty to see how a long re-interacting phase 
would produce such a result: Resonances whose interaction cross-section is small w.r.t. 
the timescale of collective expansion would generally be depleted by the dominance 
of rescattering over regeneration processes at the detectable (on-shell) mass range. 
More strongly interacting resonances would be re-thermalized at a smaller, close to 
thermal freeze-out temperature. Both of these scenarios would generally result in a 
suppression, rather than an enhancement, of directly detectable resonances. Transport 
model studies done on resonances generally confirm this [281 1201 

The observation of an enhanced ^"""/i" continuum around the p peak pUj has been 
pointed to as evidence of p broadening, which in turn would signify a long hadronic 
re-interaction phase (HOI. The absence of a broadening in the nominal peak itself, 
prevents us from considering this as the unique interpretation of experimental data. 
Moreover, even a conclusive link of broadening with hadronic re-interactions would 
still give no indication to the length of the hadronic rescattering period. Nor it would 
resolve the discrepancies pertaining hadronic resonances encountered in the previous 
paragraph; The lack of modification in either mass or width, between p-p and Au-Au 
seen so far [201 12Z] is only compatible with the NA60 result| provided the p is very 
quick to thermalize, so [201 sees only the p s formed close to thermal freeze-out. But, 
as argued in the previous paragraph, that under-estimates the abundance of the p 
and other resonances, since the pj-K^K* jK^ and even A(1520)/A ratios point to a 
freeze-out temperature significantly above the 100 MeV, commonly assumed to be the 
"thermal freeze-out" temperature in a staged freeze-out scenario. 

It is apparent that the data surveyed above is insufficient for a conclusive 
analysis. One or more missing observables are required. It is desirable that these 

I Unless one assumes that at freeze-out the SPS system to be very different from RHIC's. But no 
soft data so far would lend support to such a hypothesis, and in all statistical models mentioned here 
conditions at SPS and RHIC are comparable, save for the larger SPS chemical potentials 
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observables have a different dependence on the non-equilibrium parameters 7^^^ (see 
Eq. EJ, system volume, and temperature, so as to decorrelate these parameters and 
falsify the non-equilibrium picture in [S( 117 and it's explanation for the "horn" and 
related discontinuities. It is also desirable that these observables be strongly ensemble 
Dependant, to test the Canonical origin theory of strangeness enhancement. Finally, 
a direct experimental sensitivity to the re-interaction period would be invaluable to 
resolve the conundrums of the preceding paragraphs. 

As we will show, considering both yields and fluctuations within the same analysis 
could have all three requirements. 



2. Fluctuations in the statistical model 



Event-by-event particle fluctuations have elicited theoretical |^ |321 EHl IM] and 
experimental interest E^l EH EHI i both as a consistency check for existing models, 
as a test of hadron gas equilibration and as a way to search for new physics. The ability 
of the SHM to describe not just averages, but event-by-event multiplicity fluctuations 
has however not been widely investigated, and its applicability is currently a matter 
of controversy. 

In this work, we use the Grand Canonical ensemble (most physically appropriate 
for describing RHIC data implemented in a publicly-available software package 

1211 EHI to calculate fluctuations ((AiV)^) and yields {{N)). At freeze-out, these are 
given by 

(^.> ^ J-^^^ (1) 

((AA^.)'>-A.^(iV,) (2) 

The parameter Ai corresponds to the particle fugacity, related to the chemical potential 
by \i = e^^l^ . Provided the law of mass action holds, it should be given by the product 
of charge fugacities (flavor, isospin etc.). It is then convenient to parametrize it in 
terms of equilibrium fugacities A"'' and phase space occupancies 74. For a hadron with 
g(g) light quarks, s(s) strange quarks and isospin J3 the fugacity is then 

A.-Ar7r'7:+^ , \T^y\-~'K'~'A\ (3) 

If the system is in chemical equilibrium then detailed balance requires that 7g = 7s = 
1. Assuming 7^ = 1 and fitting particle ratios gives the T ~ 160 — 170 MeV freeze- 
out temperature typical of chemical equilibrium freeze-out models at SPS and RHIC 

ElEllHI. 

In a system expanding and undergoing a phase transition, however, the condition 
of chemical equilibrium no longer automatically holds, so one has to allow for the 
possibility that 7^ ^ 1, 7^ 7^ 1. In particular, if the expanding system undergoes a fast 
phase transition from a QGP to a hadron gas, chemical non-equilibrium ^ and super- 
cooling |41| can arise due to entropy conservation: By dropping the hadronization 
temperature to ~ 140 MeV and oversaturating the hadronic phase space above 
equilibrium (7^ ~ 1.5,7s 2), it is possible to match the entropy of a hadron 
gas with that of a system of nearly massless partons In such a case, the pion 
chemical potential is very close to, but not quite reaching the point of Bose-Einstein 
condensation Ai e™'/"^. 
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The onset of these conditions (and of the associated strangeness over-saturation 
condition 7^ > 1) is the basis of the non-equiUbrium statistical model description of 
the "horn" mentioned in the preceding section [S]. 

This limit should result in enhanced fluctuations |34| . since higher order terms 
in the Bose-Einstein expansion are always more significant than for yields: for 
AiC™'/^ = 1 — e {Ni) converges but (ANf) diverges as ~ e^^/^. Thus, a 
measure of both fluctuations and yields could serve to falsify the non-equilibrium and 
equilibrium statistical model scenarios: T and 7^ correlate in yields and anti-correlate 
in fluctuations. Measuring both yields and fluctuations can be used to determine both 
temperature and 7^ precisely. 

The final state yield of particle i is computed by adding the direct yield and all 
resonance decay feed-downs. 

(iV.)total ={m+ BjMN,) (4) 

all j—>-i 

{(AN^f) = B,^,{l-B,^,){N,)+ Y Bj^,{{AN,f) (5) 

all j — -'i all j — >i 

assuming no primordial correlation between i and j. 

The fluctuation of a ratio N1/N2 can be computed from the fluctuation of the 
denominator and the numerator [3T] (c^ = {{AX)'^'^ / {X)): 

_ ((AiVi)^) , {{AN^f) (AiViATVa) 

This is a a useful experimental observable, since volume fluctuations, difficult to model 
from first principles, cancel out event-by-event. 

Note that fluctuations of ratios scale as the inverse of the absolute normalization, 
'^Ni/N2 ~ ((^) T^)^^- Hence, a fit including such a fluctuation needs to have absolute 
normalization as one of the flt parameters. Fitting a particle multiplicity, ~ {(V) T^) 
together with the ratio provides a way to measure the system volume at chemical 
freeze-out with great precision, due to the very different scaling of the two types of 
observables. 

Note also the appearance of a negative correlation term between A^i and N2 
stemming from a common resonance feed-down (A — > pn will be a source of correlation 
between Np and A^^). 

(A7ViAiV2)«^ 6,^12 (iV,) (7) 
j 

This correlation term is an invaluable phenomenological resource, since it is sensitive 
to resonance abundance at chemical freeze-out |31) : Further reinteraction, provided 
resonances are not kicked out of the detector's acceptance region [31], preserve the 
correlation even when the original resonance stops being reconstructible. Thus, a 
measurement of the correlation together with the directly detectable resonance is a 
direct measurement of the amount of re-interaction between the thermal and chemical 
freeze-out. 

As can be seen, fluctuations measurements fulfill all three of the requirements 
enumerated in the last section. They are sensitive to the degree of chemical non- 
equilibrium and the system volume. They can gauge the amount of hadronic re- 
interaction after hadronization. And, as shown in j42j . they are strongly ensemble- 
specific. In the next section, we will attempt to interpret current RHIC data in both 
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Figure 1. Left panel: Statistical significance profile for the 200 GeV fit including 
(solid line) and excluding (long-dashed line) the cr"^"^ data-point. Right panel: 
Parameter sensitivities for the fits for T,d {V) /dy and 7s 



yields and fluctuations in light of the controversies described in the previous section, 
and to suggest experimental observables capable of further clarifying the situation. 

3. Statistical model comparison with preliminary RHIC data 

The main experimental problem with fluctuation measurements is the vulnerability 
to effects resulting from limited detector acceptance. This difficulty can be lessened, 
to some extent, by considering "dynamical" fluctuations, obtained by subtracting a 
"static" contribution which should be purely Poisson in an ideal detector. 



a'^''^^<j'~al,, (8) 

o'stat, usually obtained through a Mixed event approach |48j . includes a baseline 
Poisson component, which for a ratio N1/N2 can be modeled as 

2 _ 1 1 

'^stat 'J^ + J^ 

as well as a contribution from detector efficiency and kinematic cuts. Provided 
certain assumptions for the detector response function hold (see appendix A of JlBJ), 
subtracting astat from a should yield a "robust" detector-independent observable. 

We have performed a fit incorporating all ratios from [^, with the exception of 
the A++/p, which the STAR collaboration has taken back in more recent analyzes 
| 44) . In addition, we have included the preliminary value for the K/ir event- by-event 
fluctuation measured by STAR [37], as well as the pubhshed yield for (j) gS], tt~ 
K* (2^1 , A(1520)§ and S* [27|. As was pointed out in gT], the statistical model can 
not at present describe the difference between and cr^"^^+, although each 

§ This was taken at 5% rather than 10% centrality. The systematic error on normalization should 
not however be a dominant effect 
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Figure 2. Left panel: Fit of preliminary 200 GeV data, including the K^/tt^ 
fluctuations and the K*{&92) and A* (1520) resonance. Black circles represent a 
flt where 7q was fitted. Red triangles down, 7q = 1. Blue squares, 7q = 1 and 
fluctuation not fitted. Right panel: Predictions for resonances and fluctuations 
calculated with best fit values. Filled circles represent full resonance correlation. 
Hollow circles assume 20 % of correlations due to resonances survive within the 
detector acceptance region, compatible with the calculation in 1341 for Ar} = ±0.5 



(,s-5)=0±0.01 M 



= 0.4 ±0.01 

Au 



of the data-points can be fitted acceptably on it's own, with a parameter variation 
of 10% or so. Since these data points do not at present have a systematic error 
estimate, we do not regard this as evidence against our model. We choose, therefore, 
to fit the CT^"/^- ■ Our fit parameters include the normalization (hopefully related 
to the system "volume" at chemical freeze-out), temperature, \q,s,ii and 7q,s- We 
also require, by implementing them as data-points, strangeness, charge and baryon 
number conservation 

The results of the analysis are shown in Figs and [21 Fig. ^ (left panel) shows a 
statistical significance profile in 7^ . At each point in the abscissa, a fit was performed 
minimizing all the other parameters (T, Ag.sja, 7s, d {V) /dy). The ordinate shows the 
statistical significance of such a best fit. The right panel shows what the best fit values 
of each parameter are at each point in 7,. The left panel of Fig. |2] shows what the 
best fit to data looks like in three cases: 7, fitted (circles), 7, = 1 (triangles), and 
7q = 1 and ct^/tt not fitted (hollow squares). 

As can be seen in both Fig. nand[21(left panels), incorporating fluctuations, yields 
and ratios rules out equilibrium and forces 7^ to be well above unity. An equilibrium 
fit of the particle yields and ratios gives an acceptable description of particle ratios, 
somewhat over-predicts the K* and A(1520), but under-predicts cr^^'^ by many 
standard deviations. Choosing the Canonical ensemble for strangeness would only 
make this disagreement worse ^Sj- If ^^/^ is fitted together with ratios (red triangles 
down in Fig. |2Jl, it forces the system volume d{V) /dy to be unrealistically small 
500/to'^), thereby under-predicting particle yields by several standard deviations. 
It is only the addition of jq (circles in Fig [21) that allows fluctuations to be driven to a 
high enough value while maintaining sufficiently high volume to describe the particle 
multiplicities, and sufficiently high temperature to describe ratios. 
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The correlation between T and 7^ given by particle ratios then fixes temperature 
to T = 140 MeV, as predicted in a super-cooling scenario [SJ [2] and d{V) jdy to 
^ 1000 — 1500/m^, in agreement with earlier estimates TU' and combatible with HBT 
measurements |20) . 7s is also required to be significantly above 1, also in agreement 
with earlier studies |U1 and with a scenario where strangeness enhancement occurs 
because of the appearance of a phase where equilibrium strangeness content is higher 
than for an equilibrium hadron gas ^5 . Quantitative disagreements w.r.t. P?l ll()| in 
diy) /dy,^s^ and chemical potentials arise from data choice, and can therefore be 
considered a systematic error. 

It is important to underline that both yields and fiuctuations contribute to such 
a precise determination: Equilibrium statistical models can describe most yields and 
ratios acceptably with 7^ = 7s = 1, but fail to describe the event-by-event fluctuation. 
Conversely, transport models provide an acceptable description of event-by-event 
fluctuations j but fail to describe the yield of multi-strange particles |5U| . 

Unlike what is sometimes asserted, the A(1520) and K* are acceptably described 
by the statistical model. The strongest disagreement arises from E* under-prediction, 
at the level of 1.5 standard deviations. We await for more measurements of resonances 
such as A and p in central collisions before trying to interpret this under-prediction. 

The acceptable description of the A(1520) and K* yield using the same freeze-out 
temperature as the stable particles, and the under-prediction of the E*, makes a case 
for the proposition that the re-interaction period between hadronization and freeze- 
out might be not as signiflcant as generally thought. However, the current data is not 
capable to rule out such a long-reinteraction period. 

The forthcoming measurement of ratios such as cr'^-/^+ and f^"/^- , however, 

will make such a falsification possible, since, unlike in the case, a sizable 

resonance correlation term exists. For example, comparing 
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with the directly measured gives the amount of re-interaction Ks and tts went 
through between hadronization and freeze-out. The biggest systematic uncertainty of 
this measurement is likely to come from the need to estimate aK'-tK-Tr, the probability 
that the decay products of will be in the detector acceptance region provided the 
original resonance was. In we discuss this number, and give some estimates of it's 
effect. A more precise estimate can be obtained via the Monte-Carlo based algorithms 
used to estimate the experimental acceptance as a function of phase space for the 
directly detected K*. We believe that, while experimental acceptance has to be taken 
into account, it's effect is not strong enough to destroy the possibility of an estimate 
of 



J chemicaL j.o. 

Fig. 121 (right panel) develops further the reasoning above, and uses the best fit 
parameters obtained using the current fit to calculate both resonances and fluctuations 
likely to be measured soon. The fluctuation calculation is given with both all 
aAB^c = 1 for all A, B, C (filled circles) and qas^c = 0.2 for all A, B, C (hollow 
circles) to give an estimate of the systematics arising from resonance decays within a 
small experimental acceptance. 
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Perhaps the measurement most likely to constrain the 7g > 1 scenario is , 
since this data-point is completely independent of pionic observables. As explained in 
the preceding section, 7^ > 1 is better able to describe fluctuations because the high 
pion chemical potential gives pion fluctuations an additional boost. Thus, fluctuations 
of quantities not involving pions should be considerably suppressed since, even at 
7s ~ 2, B-E contributions for Kaons are too small to lift Kaon fluctuations too much 
above the Poisson limit. This explains why the lowest value for a'^v'^ in the given 
sample is cr^+//^-, even through resonance correlations for this ratio are not nearly 

as numerous as for the case of, for example CpJ^- It will shortly be clear weather the 

smaller (t^"^^_ is experimentally confirmed. 

In conclusion, we have used preliminary experimental data to show that, at 200 
GeV, the equilibrium model is unable to describe both yields and fluctuations within 
the same statistical parameters. The non-equilibrium model, in contrast, succeeds in 
describing almost all of the yields and fluctuations measured so far at RHIC, with the 
parameters expected from a scenario where non-equilibrium arises through a phase 
transition from a high entropy state, with super-cooling and oversaturation of phase 
space. We await more published data to determine weather the non-equilibrium model 
is really capable of accounting for both yields and fluctuations in all light and strange 
hadrons produced in heavy ion collisions. 
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Tomlinson foundation. We would like to thank J. Rafelski, C. Gale and S.Jeon for 
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